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ABSTRACT 

An impact-type Joss-Waldvogel disdrometer (JWD), a two-dimensional video disdrometer (2DVD), and 
a laser optical OTT Particle Size and Velocity (PARSIVEL) disdrometer (PD) were used to measure the 
raindrop size distribution (DSD) over a 6-month period in Huntsville, Alabama. Comparisons indicate event 
rain totals for all three disdrometers that were in reasonable agreement with a reference rain gauge. In 
a relative sense, hourly composite DSDs revealed that the JWD was more sensitive to small drops (<1 mm), 
while the PD appeared to severely underestimate small drops less than 0.76 mm in diameter. The JWD and 
2DVD measured comparable number concentrations of midsize drops (1-3 mm) and large drops (3-5 mm), 
while the PD tended to measure relatively higher drop concentrations at sizes larger than 2.44 mm in di- 
ameter. This concentration disparity tended to occur when hourly rain rates and drop counts exceeded 
2.5 mm IV 1 and 400 min , respectively. Based on interactions with the PD manufacturer, the partially in- 
homogeneous laser beam is considered the cause of the PD drop count overestimation. PD drop fall speeds 
followed the expected terminal fall speed relationship quite well, while the 2DVD occasionally measured 
slower drops for diameters larger than 2.4 mm, coinciding with events where wind speeds were greater than 
4 m s _1 . The underestimation of small drops by the PD had a pronounced effect on the intercept and shape of 
parameters of gamma-fitted DSDs, while the overestimation of midsize and larger drops resulted in higher 
mean values for PD integral rain parameters. 


1. Introduction 

Knowledge of the raindrop size distribution (DSD) is 
essential for many applications in various disciplines of 
Earth sciences. For example, rainfall is a key compo- 
nent of Earth’s water and energy cycles and its intensity 
and duration are of interest to a wide variety of Earth 
science and engineering disciplines. While rain gauges 
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can directly measure rain intensity and accumulation 
at a point, regional to global mapping of rainfall relies 
heavily on measurements that include ground-based 
radar and satellite remote sensing. These remote sensing 
measurements necessarily employ empirical relations 
and underlying assumptions that are directly related to 
and affected by characteristics of the DSD. 

The accuracy of the rainfall estimates using single- 
polarization radar measurements traditionally relies 
on the appropriateness of derived radar rainfall rela- 
tions including radar reflectivity-rain-rate relationships 
(i.e., the Z-R relation). Both radar measurements and 
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rainfall are integral products of the DSD and a surface- 
based disdrometer is often employed to derive Z-R 
relations for a climate region, a particular weather system 
(e.g., frontal and tropical), or even a segment of a storm 
(Atlas et al. 1999; Radhakrishna and Narayana Rao 
2010). The highly variable nature of the DSD, for ex- 
ample from one weather system to another or between 
convective and stratiform segments of the same storm, 
results in substantial differences in derived Z-R relations. 
Considering the disdrometer-based Z-R relations, the 
variability in the Z-R also depends on the choice of dis- 
drometers (e.g., impact type and optical), method of 
derivation (e.g., linear versus nonlinear least squares fit), 
and the data processing techniques (e.g., time versus rain 
rate-based averaging DSD observations) (Campos and 
Zawadzki 2000; Tokay et al. 2001; Lee and Zawadzki 
2005a, b). 

For the dual-polarization radars, disdrometer obser- 
vations are often employed in deriving relationships 
between polarimetric radar observables such as reflec- 
tivity at horizontal polarization (Zh), differential reflec- 
tivity (Zdr), and specific differential phase (Kdp) and rain 
rate (e.g., Bringi and Chandrasekar 2001; Tokay et al. 
2002; Brandes et al. 2002; Bringi et al. 2003; and many 
others). Based on combinations of the polarimetric vari- 
ables, a hybrid algorithm can then be developed for es- 
timating rainfall (Cifelli et al. 2002; Bringi et al. 2002). 
The quality of polarimetric radar observables is critical 
for the accuracy of rainfall estimation. One of the nu- 
merous methods to quality control radar measurements 
uses a disdrometer-derived consistency equation among 
the three variables (Scarchilli et al. 1996; Ryzhkov et al. 
2005). 

Disdrometer observations are often employed to de- 
rive the parameters of model DSDs. The parametric 
form of the DSD (e.g., gamma, exponential, or lognor- 
mal distribution models) plays an important role in both 
active and passive satellite-based microwave sensor 
rainfall estimation. For example, the National Aeronautics 
and Space Administration (NASA)’s Global Precipita- 
tion Measurement (GPM) mission will launch its core 
satellite in 2014, carrying a dual-frequency precipitation 
radar (DPR) (Hou et al. 2008). In fact, efforts are al- 
ready underway to demonstrate the DPR algorithm for 
the dual-frequency retrievals of the DSD (Mardiana 
et al. 2004; Meneghini and Liao 2007). In a related fash- 
ion, there have been multiparameter radar-based efforts 
to determine DSD variables from dual-polarization 
radar (Bringi et al. 2003, 2002; Zhang et al. 2001; Thurai 
et al. 2012) and from vertically pointing radars (Williams 
et al. 2000; Cifelli et al. 2002). The main motivation in 
these studies is to better estimate the surface rainfall 
using a parametric form of the DSD and more accurate 


estimation of the parameters associated with the par- 
ticular DSD model used. 

Disdrometer measurements are often used to de- 
termine the accuracy of the DSD retrieval from space- 
borne radar measurements. Munchak and Tokay (2008), 
for instance, employed the impact-type Joss-Waldvogel 
disdrometer (JWD) observations to simulate the re- 
trieval of parameters of the DSD from DPR. They used 
a three-parameter gamma distribution to represent the 
DSD in mathematical form. The gamma distribution has 
also been employed by the NASA’s Tropical Rainfall 
Measuring Mission (TRMM) precipitation radar (PR; 
Iguchi et al. 2000; Kozu et al. 2009a, b) and the GPM 
DPR (Nakamura and Iguchi 2007) algorithm teams. The 
gamma-model DSD is expressed as 

N(D) = N 0 D' n exp(— AD), (1) 

where N 0 , m, and A are intercept, shape, and slope pa- 
rameters, respectively and D is the drop diameter. Kozu 
et al. (2009a) employed Z = ARP relations, where A is 
the coefficient and b is the exponent, which were origi- 
nally derived mostly from disdrometer observations 
in the tropics and produced a global average Z-R for 
convective and stratiform rain. They then derived the 
power-law relation between slope and intercept pa- 
rameters for convective and stratiform rain by fixing 
m = 3. There is considerable interest in eliminating as- 
sumptions of a constant shape parameter in ground- 
based radar and GPM DPR algorithms alike. In this 
regard, considerable effort has been made in the past 
and present to find robust relationships between the 
three parameters of the gamma distribution (Zhang et al. 
2001; Brandes et al. 2002; Vivekanandan et al. 2004). To 
this end, long-term observations of the disdrometer ob- 
servations in various climate regimes are useful to test the 
robustness of the relations between the DSD parameters. 

Disdrometer observations are useful for studying the 
highly variable nature of the DSD, at least at the point of 
measurement. Disdrometers can provide a temporally 
continuous record of the DSD and can address temporal 
variability of rainfall between consequent radar scans 
and satellite overpasses. The satellite- and radar-based 
DSD retrievals implicitly assume the homogeneity within 
the footprint and pixel space, respectively. Tokay and 
Bashor (2010) employed three JWD observations to 
study the variability of DSD and integral rain param- 
eters within the radar pixel. The efforts are currently 
underway to study the DSD variability within TRMM 
PR footprint using JWD and laser optical Particle Size 
and Velocity (PARSIVEL) disdrometer (PD). 

Disdrometer observations have a wide range of ap- 
plications and hence users from various disciplines have 
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a vested interest in understanding measurement accu- 
racy and instruments’ shortcomings. While accurate 
measurements of raindrop size and number character- 
istics across the full range of drop sizes are ideal, a par- 
ticular drop size range may have higher importance for 
a specific application. In fact, the design and develop- 
ment of a disdrometer is sometimes driven by a partic- 
ular measurement need. The JWD, for instance, was 
developed to measure radar reflectivity and therefore 
the measurement accuracy of midsize (1-3 mm in di- 
ameter) and large (larger than 3 mm in diameter) drops 
was relatively more important than measuring small 
drops (less than 1 mm in diameter; Joss and Waldvogel 
1967). As another example, the ability to image in- 
dividual hydrometeors shape and fall velocity were 
a primary motivation for the development of the two- 
dimensional video disdrometer (2DVD; Schonhuber 
et al. 2007). The identification of the precipitation type, 
measurement of size, and fall velocity of the hydrome- 
teors, while keeping the costs low, were the main themes 
behind the development of the PD (Loffler-Mang and 
Joss 2000). 

The three disdrometers (JWD, 2DVD, and PD) used 
for the present study are commercially available and 
have been used extensively in the literature. Neverthe- 
less, there are other types of optical disdrometers that 
are either commercially available or are just developed 
for scientific use. Low-power radars can also be used as 
radar disdrometers even though they do not directly 
measure the individual hydrometeors. One example, the 
Precipitation Occurrence Sensor System (POSS), is a 
bistatic continuous wave X-band radar that was origi- 
nally developed as an automated present weather sensor 
(Sheppard 1990). However, it has been widely used to 
measure the DSD (e.g., Lee and Zawadzki 2005a,b) 
and its small-scale spatial variability (Lee et al. 2009). 
Miriovsky et al. (2004) intended to quantify the spatial 
variability of radar reflectivity within the approximate 
spatial scale of a radar beam (~1 km 2 ) employing POSS, 
2DVD, JWD, and two optical spectropluviometers 
(Hauser et al. 1984). They concluded that the inter- 
instrument differences/errors in measuring the capa- 
bility of the DSD did not allow determination of 
small-scale variability of radar reflectivity. 

There have been a number of field studies designed to 
evaluate the performance of a disdrometer. Since each 
disdrometer has its own shortcomings, collocated ref- 
erence rain gauges are often collocated with the dis- 
drometers. Outside of the systematic errors mainly 
because of the wind, collocated gauges provide reliable 
rain totals (Tokay et al. 2010). However, good agree- 
ment between gauge and disdrometer rain totals does 
not guarantee good performance of the disdrometer. 


For example, small drops often do not substantially 
contribute to the rainfall and therefore the instrument 
errors in these size ranges are masked in comparisons. 

During TRMM field campaigns in east-central Flor- 
ida and the Amazon basin of Brazil, the JWD and 2DVD 
were collocated to obtain higher accuracy and to ascer- 
tain uncertainty in the DSD estimates for radar rainfall 
estimation. While the disdrometers exhibited reasonable 
agreement with collocated rain gauge totals, there were 
noticeable differences between the JWD and 2DVD in 
the small and large size ends of the drop size spectrum 
(Tokay et al. 2001, 2002). Krajewski et al. (2006) conducted 
a field campaign where 2DVD, PD, and a dual-beam 
spectropluviometer (DBS) were collocated. Similarly, 
the disagreement between the disdrometers was at both 
ends of the drop size spectrum. While the differences in 
small drops were attributed to the measurement errors, 
differences in large drop counts were attributed to sam- 
pling errors. 

Optical disdrometers can also measure the fall velocity 
of raindrops. Krajewski et al. (2006) reported that the 
majority of the fall speed measurements followed Gunn 
and Kinzer (1949, GK hereafter) terminal fall speed in 
DBS, but drops larger than 2 mm in diameter de- 
viated from GK terminal fall speed in 2DVD and PD 
measurements. 

The highly variable nature of the DSD results in 
sampling errors in disdrometer measurements. Time 
averaging can reduce the sampling errors, but it may also 
mask the intrinsic physical variability of the DSD. Col- 
locating several of the same type of disdrometers as in 
the TRMM field campaigns can help to determine the 
sampling error even though rain is not a homogeneous 
process and each instrument has its own calibration. As 
such, Tokay et al. (2005) showed noticeable differences 
in rainfall estimations, which were derived from six col- 
located JWD observation-based radar rainfall relations. 
Cao et al. (2008), on the other hand, employed two 
collocated 2DVD observations to refine the relation be- 
tween the slope and shape parameters of the gamma 
distribution. Reducing the sampling errors, they ob- 
tained a good agreement between radar and disdrometer 
DSD and rain parameters. Thurai et al. (2011) examined 
several rain events where two 2DVDs and two PDs were 
collocated. While there was a good agreement in derived 
DSD and rain parameters between the same types of 
disdrometers, PD had significantly higher mean mass 
diameter and rain rate than 2DVD, particularly at rain 
rates above 30 mm 1C 1 . 

In the last decade, there have been numerous de- 
velopments in commercially available disdrometers, 
ostensibly to improve the breadth of measurement accu- 
racy. Accordingly, this study provides a new comparison 
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of the disdrometers’ performance. This type of compar- 
ison is important because a complete understanding of 
instrument measurement uncertainty provides a frame- 
work for assessing the uncertainty in any retrieval of 
DSD, against which the disdrometers information is 
compared. Specifically, 6 months of collocated obser- 
vations from JWD, 2DVD, and PD were collected in 
Huntsville, Alabama, and have been employed to in- 
vestigate and characterize disdrometer performances 
for DSD and selected rain parameters. This study 
provides a brief description of each disdrometer type in 
section 2. Sections 3 and 4 describe the rainfall statistics 
and the observation site and data, respectively. A com- 
parison of event rain totals between the disdrometers 
and the collocated rain gauges is presented in section 5, 
while hourly DSD characteristics observed by the dis- 
drometers are given in section 6. Section 7 shows a com- 
parison of fall velocity measurements from the 2DVD 
and PD. A comparison of the selected DSD and rain 
parameters is demonstrated in section 8 and conclusions 
are presented in the last section. 

2. Disdrometers 

a. ./ oss-Wald vogel disdrometer 

The impact-type JWD consists of sensor head and 
signal processing electronics (Joss and Waldvogel 1967). 
The sensor head has a cone, which is a Styrofoam body 
covered with an aluminum foil on the hitting surface. 
The cone has been equipped with wires, which are put 
on the hitting surface to help water drops leaving the 
surface and thus avoiding the formation of puddles. 
The RD-80 is a signal processing box that integrates the 
analog signal processing circuits and analog-to-digital 
(AD) converter into one unit. Prior to 2000, the processor 
box (RD-69) had an analog output, which was connected 
to the AD (M. Weibel, Distromet, Ltd., 2011, personal 
communication). To position the sensor away from taller 
structures, a 100-m detachable cable can now be used to 
connect the sensor with the indoor signal processing 
box. 

JWD measures the raindrop size from 0.3 to 5.5 mm in 
diameter with an accuracy of ±5% if the drops are 
evenly distributed over the sensitive surface. The stan- 
dard output includes drop counts at 20 uneven channels 
where the channel width ranges from 0.1 to 0.5 mm at 
1-min intervals. A more detailed output with 127 size 
channel drop counts at 10-s intervals is also available. 
The manufacturer provides a calibration table where the 
accuracy of 127 size bins is presented. The Distromet, 
Ltd. company has been manufacturing the JWD since 
1972 and, based on our knowledge, the JWD has been 


used more than any other disdrometer in the literature. 
Additional information about the JWD can be found at 
the company website (http://www.distromet.com). 

b. Two-dimensional video disdrometer 

The 2DVD consists of an outdoor sensor unit contain- 
ing dual-optical line scan cameras and an indoor user 
terminal. The current third-generation 2DVD, which is 
also referred to as the compact version because of its 
smaller dimensions and fewer number of components 
than previous versions, was first delivered to the GPM 
ground validation group in 2009. We are one of the first 
to use data from this compact 2DVD in addition to 
a companion study reported in Thurai et al. (2011). 
Note that the compact 2DVD was designed as an im- 
provement to the second-generation low-profile unit 
(single optical bench, more rugged instrument) and 
that both the compact and low-profile units were 
developed to improve on the original taller 2DVD 
(Kruger and Krajewski 2002). Here, the compact and 
low-profile designs were driven by the need to reduce 
wind-induced measurement errors in response to Nespor 
et al. (2000). Like JWD, the 2DVD requires shelter and 
power. 

The 2DVD measures the hydrometeor size, fall ve- 
locity, and shape through its two high-speed line scan 
cameras. The sensing cross section is approximately 
100 cm 2 . The imaging grid resolution for raindrops is less 
than 0.2 mm in diameter in both horizontal and vertical 
directions. Each hydrometeor is measured twice through 
two orthogonally oriented cameras situated in offset 
measurement planes that are separated by approximately 
6 mm in a vertical direction. The optical alignment and 
the hydrometeor matching are critical for the accuracy of 
the measurements. The 2DVD is calibrated by dropping 
calibration spheres ranging from 0.5 to 10 mm in di- 
ameter. The raw output includes the time stamp, equiv- 
alent diameter, measured fall velocity, oblateness, and 
height, width, and position of the hydrometeor at each 
plane. Among other disdrometers, the 2DVD provides 
the most detailed information about the individual hy- 
drometeors and the 2DVD data viewer menu presents 
the image of the hydrometeors in two measuring planes. 
Schonhuber et al. (2007) summarized all the technical 
improvements of the 2DVD until its publication and 
more detailed information can be found on the manu- 
facturer’s (Joanneum Research in Graz, Austria) website 
( www. distrometer .at) . 

To improve data quality beyond what the manufac- 
turer’s software provides, drops exceeding ±50% of 
their theoretical terminal fall speed (Beard 1976) were 
considered either secondary (i.e., a result of splash) or 
mismatched drops and eliminated from the dataset, 
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which resulted in an elimination of 8% of the observa- 
tions. Whereas observations obtained with the first- 
generation 2DVD during TRMM field campaigns in 
central Florida and the Amazon basin of Brazil, 18% 
and 22% of the observations, respectively, were elimi- 
nated after applying the same quality-control method. 
Thus the faster cameras and new design reduced the 
production of secondary and mismatched drops. 

c. PARSIVEL disdrometer 

The PD is a laser optical disdrometer and was con- 
structed by PM Tech AG, Pfinztal, Germany. The PM 
Tech sold the right of the PD to OTT Hydromet, Kempten, 
Germany, after four years of operation in 2004 
(M. Loffler-Mang 2007, personal communication). OTT 
Hydromet improved the calibration of the disdrometer 
but used a more inexpensive laser device to reduce the 
cost substantially. Indeed, since it started its production, 
OTT Hydromet has sold approximately 1000 PD units 
(K. Nemeth, OTT, 2011, personal communication), many 
of which are used as present weather sensors. 

The PD measures the size and fall velocity of the 
hydrometeors (Loffler-Mang and Joss 2000). The raw 
output is the number of drops in 32 size and 32 fall ve- 
locity categories. The size and fall velocity ranges are 
from 0.2 to 25 mm and from 0.2 to 20 m s _1 , respectively. 
The first two size categories, which correspond to sizes 
less than 0.2 mm, have been left empty because of the 
low signal-to-noise ratio. The OTT PD laser beam is 
180 mm long, 30 mm wide, and 1 mm thick. The nominal 
cross-sectional area is 54 cm 2 , but the effective sampling 
cross section including the edge effects is expressed as 
180 X (30 — L/2), where L is the size parameter. The 
raindrops are approximately spherical for sizes less than 
1 mm in diameter, and therefore the size parameter is 
the equivalent diameter for raindrops below this size. 
For larger raindrops, a correction for oblateness has 
been made to the measured size parameter. The PD 
measures the fall velocity by monitoring the duration of 
the hydrometeor presence in its laser beam. A correc- 
tion has been applied to the fall velocity for size pa- 
rameters between 1 and 5 mm (M. Loffler-Mang and 
U. Blahak 2005, personal communication). Battaglia 
et al. (2010) presented the measurement principals and 
accuracies of the PD size and fall velocity measure- 
ments for snowflakes. The PD data output is provided 
either at 10-s or 1-min intervals. 

In addition to the data processing discussed above, 
which was conducted internally by the manufacturer’s 
software, and to be consistent with the 2DVD observa- 
tions, we removed drops from our PD dataset that ex- 
ceeded ±50% of their theoretical terminal fall speed 
(Beard 1976). This criterion eliminated 20% of the drops. 


3. Rainfall statistics 


The percent bias and absolute percent bias are used 
to evaluate the performance of rain gauges and dis- 
drometers. The percent bias and absolute bias between 
the two instrument measurements ( x , y) for n samples 
are calculated as 

, • bias 

percent bias = -, and (2) 

\x,y) 


, , , . absolute bias 

percent absolute bias = r 

{x,y) 


(3) 


where bias, absolute bias, and the mean value between 
the two variables (x, y) are expressed as follows: 

1 11 

bias = -£(*,-%■), (4) 

n (= j 


absolute bias = - lx. — y\, and 
n i=1 

1 A 

( x >y)=z X 


i=t 

_i v (*«■ + yJ 


i=t 


(5) 

( 6 ) 


4. Measurement site and data 

A rain-measuring instrument test site was constructed 
in Huntsville, Alabama (34.72°N, 86.64°W). This study 
employs NASA’s newly acquired third generation of the 
2DVD, a JWD, and an OTT PD (Fig. 1). The field study 
was between 18 December 2009 and 2 June 2010. 

A Texas Electronics (model TR-525I) tipping-bucket 
(TB) gauge was operated throughout the experiment 
and three MetOne Inc. tipping-bucket gauges were 
deployed near the end of the field study. One of the 
MetOne gauges malfunctioned as a result of low battery 
in the datalogger system while the other two MetOne 
gauges exhibited excellent agreement with 0.4% bias 
and 1 % absolute bias for the 22 rain events (Fig. 2a). The 
agreement between the two MetOne gauges and the 
Texas Electronics gauges was very good with 5 % bias and 
8% absolute bias between the gauge types (Figs. 2b,c). 
Because it was located at the site the longest, the Texas 
Electronics gauge was considered the reference gauge for 
evaluating the disdrometers performance. The site was 
also equipped with an anemometer situated on a 10-m 
tower, which allowed us to investigate the wind effect on 
size and fall velocity of raindrops. 


5. Rainfall measurements 

Throughout the field study, the 2D VD recorded 374 mm 
of rainfall in 7054 rainy minutes over 32 rain events. A 
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Fig. 1. In situ precipitation measuring test site. The instruments 
used in this study are labeled as (top) JWD sensor inside the white 
painted box, TB rain gauge back to the right and one of the sec- 
ondary TB gauges back to the left; (top), (middle) 2DVD in the 
back, and PD on the left and in front, respectively, (bottom) A 
close look at the third generation of 2DVD. 

rain event is defined here as having at least 1 mm of rain 
accumulation and individual events were segmented by 
1 h or more rain-free periods. A rainy minute must have 
at least 10 drops and a rain rate of 0.1 mmh' 1 or higher, 


similar to the quality control performed by Tokay et al. 
(2005) on JWD datasets. Integrated over all events, the 
2DVD recorded only 0.4% more rainfall than the rain 
gauge, which recorded 372 mm of rainfall. The JWD and 
PD, on the other hand, had 3.5% less and 7.5% more 
rainfall than the gauge, respectively. The JWD and PD 
also had 143 more and 117 less rainy minutes than the 
2DVD, respectively. This resulted in higher conditional 
rain rate (3.5 mmh -1 ) for the PD and lower conditional 
rain rate (3.0 mm h -1 ) for the JWD. The conditional rain 
rate is calculated here as the ratio of total accumulation 
to the rainy minutes. 

The JWD accumulated more rainfall than the gauge 
when the gauge event rain total was less than 6 mm and 
vice versa for higher rain totals (Fig. 3a). Higher rain 
total events in gauge data are usually associated with the 
occurrence of heavy rain events where the JWD's dead 
time can result in underestimation of rainfall (Tokay 
et al. 2005). An empirical dead time correction can be 
used to increase the JWD estimate of rainfall, but it was 
not used in this study since it is a multiplicative matrix 
and therefore does not add any drops when the size bin 
is empty resulting in an unrealistic size distribution. The 
2DVD also accumulated more rainfall than the gauge 
when the event rain total was less than 5 mm; however, 
there is very little disagreement between the two at 
higher rain totals (Fig. 3b). This indicates that the 
gauge may have underestimated in some events. A 
single tip of the bucket can have a more pronounced 
effect on gauge totals at light rainfall, and a partially 
full bucket that has not tipped may explain some of 
the event rain total differences. Unlike the 2DVD and 
JWD, PD had higher accumulations than the gauge in 
the majority of rain events (Fig. 3c). Among the three 
disdrometers, the 2DVD showed the best agreement 
with the rain gauge having the lowest relative and 
absolute biases. The PD, on the other hand, exhibited 
the highest biases and absolute biases relative to the 
gauge, but the absolute bias was still less than 15%, 
which has been considered an upper threshold by the 
JWD manufacturer. The JWD and 2DVD exhibited 
the best agreement among the disdrometers (Fig. 3d), 
and JWD and PD differed the most in event rain totals 
(Fig. 3e). 

6. Raindrop size distribution measurements 

As a preface to comparing specific characteristics of 
the rain DSD between the disdrometers, it is important 
to note that the agreement between disdrometer event 
rain totals suggests that their respective measurements 
of drop numbers in the midsize range are in general 
agreement; this is because rainfall is mostly contributed 
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ab. bias = 7.9% 
bias = -5.3% 



Rain Total (mm) 


Rain Total (mm) 


Rain Total (mm) 


Fig. 2. Comparison of event rain totals between two MetOne Inc. (RG1 and RG2) and the Texas Electronics (RG3) 
rain gauges for 22 rain events. The percent bias and absolute bias are also given. 


by the midsize raindrops. Now in order to more closely 
examine departures in DSD agreement among the 
instruments, we move away from rainfall accumulations 
and focus on 1-min DSD observations averaged to 1-h 
intervals. Accordingly, we have selected 70 h of rainfall 
where all three disdrometers had more than 10 rainy 
minutes of observations and no more than 3-min dif- 
ferences were observed by all disdrometers to avoid the 
sampling fluctuations. The maximum drop diameter was 
also required to exceed 3.5 mm in all disdrometer hourly 
spectra to ensure that large drops were adequately sam- 
pled. Figure 4 shows 24-hourly composite spectra from 
15 different rain events. 

Previous observations between 2D VD and J WD from 
shorter field studies showed that the disagreement be- 
tween the two disdrometer measurements occur at small 
and large drops (Tokay et al. 2001, 2002). The concave 
downward shape of JWD size spectra at the smaller drop 
end was interpreted as a physical fact as well as an in- 
strument artifact. It was considered that evaporation 
exceeds collision breakup in convective rain resulting in 
a concave downward shape of the size spectra as shown 
in western tropical Africa (Sauvageot and Lacaux 1995) 
and the western tropical Pacific Ocean (Atlas and 
Ulbrich 2000). The background noise and disdrometer 
dead time in the JWD measurements were also attrib- 
uted to the lack of small drops (Tokay et al. 2003, 2005; 
Tokay and Bashor 2010). 


In this study, the concentration increased toward 
smaller sizes, and peak concentration occurred in either 
of the first two size bins (centered at 0.35 and 0.45 mm) 
of the JWD in the majority of the hourly spectra (Fig. 4). 
At the same time, 14% of the observations had a peak 
concentration at 0.651mm, but these observations did 
not necessarily coincide with the high rain rates. This 
shows that the background noise rather than dead time 
may cause the underestimation of small drops. For 
2DVD, the peak concentration was mainly at 0.7 mm 
but also observed at 0.3, 0.5, and 0.9 mm in diameter. 
The PD had its peak concentration at 0.84 and 0.96 mm. 
This demonstrated that JWD was more sensitive to the 
small drops, and PD is the least sensitive. In a few 
spectra, the increase in concentration toward smaller 
diameters was not continuous in the JWD data (Figs. 
4t,u). The concentrations were about the same at 0.76 
and 0.9 mm for JWD. It should be noted that the JWD 
was originally developed to calculate reflectivity, which 
is proportional to the sixth moment of drop diameter 
(Joss and Waldvogel 1967). The concentration of small 
drops has a pronounced effect on total concentration 
and therefore on the modeling aspect of the DSD. 

The sampling cross section of the 2DVD is approxi- 
mately twice as large as the JWD and PD, and therefore 
it is expected that the 2DVD may sample the very large 
drops (larger than 5 mm in diameter) better than the 
other two disdrometers. The largest size bin of the JWD 
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Fig. 3. Comparison of event rain totals between rain gauge, JWD. 2DVD, and OTT PD for 32 rain events. The bias 

and absolute bias are also given. 


is centered at 5.14 mm, and drops of this size and larger 
are registered in this size bin. Since the JWD size 
measurement depends on the force at which drops im- 
pact the measurement cone and the drop fall velocities 
only incrementally increase at these very large drops, 
the JWD cannot distinguish the very large drops. The 
sampled maximum drop diameter is generally quite 
different between the disdrometers because of sampling 


limitations. Fortunately, the very large drops occur 
quite rarely, and as a result these drops typically do not 
contribute significantly to the rainfall parameters, ex- 
cept for reflectivity. Yet, the integral rain and DSD pa- 
rameters that are calculated from model spectra do 
require the knowledge of maximum drop diameter. 
The largest drop diameter ever reported in the litera- 
ture is 8 mm in diameter, which was observed below 


Concentration (m * 3 mm -1 ) Concentration (m ' 3 mm* 1 ) Concentration (m * 3 mm* 1 ) Concentration (m * 3 mm* 1 ) 
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Fig. 4. Comparison of hourly raindrop spectra from selected rain events. The hourly spectra are given for JWD (solid), 2DVD (dashed- 
dotted), and OTT PD (dashed). The hourly wind speed, total number drops min -1 , and rain rate are also given. 
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5. The difference between OTT PD and 2DVD hourly number concentration at 2.8 mm in diameter as a function 
of (a) hourly wind speed, (b) hourly total of drops min -1 , and (c) hourly rain rate. 


the cloud base using aircraft-mounted particle probes 
(Beard et al. 1986). 

The maximum drop diameter in this study was 7.74 mm 
in diameter, which was observed by the 2DVD. Although 
PD reported a number of drops at its 23rd size bin, which 
is centered at 7.72 mm in diameter, the accuracy of size 
measurements is questioned and explained later in this 
section. The PD also had the largest drop sizes in most of 
the hourly composites. The 2DVD recorded the largest 
drops in 30% of the hourly size spectra. The largest drops 
were larger than the JWD maximum drop diameter limit 
in 24% and 34% of the hourly composites in the 2DVD 
and PD, respectively. It should be noted that PD bin 
width is 1 mm at sizes between 5 and 10 mm in diameter, 
and this makes the uncertainty quite high at these sizes. 

The most striking feature among hourly composites 
was the relatively higher concentrations of drops larger 
than 2.4 mm in diameter in PD DSDs. This feature was 
visually observed in numerous composites (Figs. 4d,f,j,w) 
and cannot be attributed to inadequate sampling. Agree- 
ment between the 2DVD and JWD midsize concentra- 
tions suggest that the PD overestimates drops larger than 
2.4 mm in diameter in some cases. This feature was not 
observed when the PM Tech PD was compared with JWD 
at the NASA Wallops Flight Facility (WFF), Wallops 
Island, Virginia. This difference may be because OTT used 
a less expensive laser device, which is approximately 20% 
inhomogeneous across the beam (K. Nemeth, OTT, 2010 
personal communication). As far as an overestimation of 
large drops by the PD, Thurai et al. (2011) first reported 
this problem from a comparison of dual-2DVD and dual- 
PD mass-weighted drop diameters Z) mass at a 1-min reso- 
lution. They found that PD had a higher T> mass than the 
2DVD when £> mass was larger than 2 mm and the rain rate 
exceeded approximately 20mmh _1 . 


In this study, we used a much larger dataset from 
the same site. Rather than £> mass , the difference in 
concentration between PD and 2DVD at 2.8 mm was 
examined as a function of wind speed, drop counts, 
and rain rate. Figure 5a shows that there is no wind 
speed dependency, while Figs. 5b and 5c indicate 
a dependency on the number of drops and rain rate. If 
the concentration difference of 1 drop m -3 mm -1 and 
greater is considered to be significant at this size, the 
disdrometer observations can be categorized three 
ways based on drop counts and rain rate. If the drop 
counts are less than 400 per rainy minute, the con- 
centration difference is insignificant and vice versa 
is true at drops counts larger than 1000 per rainy 
minute. Similarly, the concentration difference is in- 
significant at rain rates less 2.5 mm h^ 1 and is signifi- 
cant at rain rates above 4.5 mmh -1 . It should be noted 
that hourly rain rates in Fig. 5c were calculated from 
only rainy minutes. Since more drops hit the sampling 
cross section during more intense rainfall, it is likely that 
a greater number of drops are subject to fall through 
the nonuniform cross section of the PD laser, thereby 
providing a greater chance for large drops to be falsely 
reported. 

7. Fall velocity measurements 

In addition to the drop size, the 2DVD and PD mea- 
sures the fall velocity of hydrometeors. After binning the 
2DVD size measurements using PD-size bins, the mean 
and standard deviation of the raindrop fall velocities were 
calculated for each size bin for 12 rain events. These fall 
velocity calculations are shown against the interpolated 
terminal fall speed of raindrops (Beard 1976) in Fig. 6. 
Each size bin was required to have at least 10 drops in 
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Fig. 6. Mean (squares) and standard deviation (vertical bars) of measured fall velocity by 2DVD as a function of diameter for 12 rain 
events. The drop counts are binned following OTT PD bin intervals. The terminal fall speed is also shown (line). 


a given event. The mean measured fall velocities match 
the terminal fall speed across all size ranges except for 
sizes above 2.44 mm in diameter, where the mean fall 
velocity was less than the terminal fall speed. The dif- 
ferences between the two fall velocities are apparent in 
several events (Figs. 6b,f-h). Figure 7 shows the differ- 
ences in fall velocities in relation to the mean event 
horizontal wind speed estimated from a collocated an- 
emometer. For a 3.35-mm drop, the terminal fall speed is 
8.4 ms -1 , and three rain events had a 6% or greater 
difference in fall speeds when the wind speed exceeded 


4 ms -1 (Fig. 7a). At 3.86 mm, the difference was 7% or 
higher (Fig. 7b). 

The PD fall velocities were higher than the inter- 
polated terminal fall speed at sizes less than 1.22 mm and 
vice versa was the case for the larger raindrops with the 
exception of very large drops, which had higher terminal 
velocities (Fig. 8). No event-to-event variability was 
observed. This may indicate that the slower fall velocities 
in 2DVD DSDs could be related to the measurement 
error. However, the measurement principles of the 
2DVD and PD are different, and the reasoning for 
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Fig. 7. Differences between terminal fall speed and measured fall velocity of 2DVD as 
a function of environment wind speed at (a) 3.35 and (b) 3.86 mm in diameter. 


2DVD fall velocity deviation from predicted terminal 
fall speed is beyond the scope of this study. It should be 
noted that the OTT PARSIVEL measures fall speeds 
more accurately than the PM Tech PARSIVEL. PM 
Tech fall velocities are higher than terminal fall speed at 
drop sizes less than 0.84 mm and are less than terminal fall 
speed at higher sizes as shown from the measurements 
during the Canadian CloudSat/Cloud-Aerosol Lidar 
and Infrared Pathfinder Satellite Observations (CALI- 
PSO ) Validation Project (C3VP) in Fig. 9. The differ- 
ence between the fall speeds exceeded 1ms -1 at 
a diameter of 4.38 mm. 

An optical spectrometer developed by Donnadieu 
(1980) and upgraded by Hauser et al. (1984) also shows 
the significantly slower fall speeds at sizes larger than 
3 mm (Salles and Creutin 2003). Salles and Creutin 
(2003) attributed the slower fall speeds to the down- 
drafts and determined the size measurements errors 
of the JWD. Krajewski et al. (2006) used an updated 
model of dual-beam spectrometer and found the 
agreement with the terminal fall speed was quite good 
for the drops less than 3 mm but there were only a few 
drops at larger sizes. They also found that the PM Tech 
PARSIVEL had slower drops at sizes above 1.4 mm in 
diameter and attributed this feature partly to the 
quantization of the velocity versus size relationship. 
The first-generation 2DVD also had slower drops ex- 
cept at small sizes in their study, and this feature was 
attributed to the wind effect around the bulky structure 
of the instrument (Nespor et al. 2000). Barthazy et al. 
(2004) developed a dual-beam hydrometeor velocity 


and shape detector and found the raindrop fall velocities 
to be in good agreement with the expected terminal fall 
speed between 1- and 2.5-mm drops. Drop splash resulted 
in lower fall speeds at the smaller drop sizes and slightly 
higher fall speeds were observed for drops larger than 
2.5 mm. 


8. Raindrop size distribution and integral rainfall 
parameters 

As described in the introduction, the three-parameter 
gamma function is often employed to parameterize the 
disdrometer-observed DSD measurements in various 
applications including precipitation retrieval algorithms 
from TRMM and GPM radar measurements. In this 
study, two different normalized gamma functions are 
employed to determine the role of different disdrometer 
measurements on DSD. The normalization was done 
with respect to the total concentration N T and liquid 
water content w (Tokay and Bashor 2010). The nor- 
malized intercept parameters with respect to total 
concentration Nf and liquid water content N w are 
expressed as 


N r 

w* = Z_ 

T D 

mass 


and 


N. = 


256 10 3 W 

17 Pw ^mass 
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Fig. 8. Mean (squares) and standard deviation (vertical bars) of measured fall velocity by OTT PD as a function of drop diameter for 12 

rain events. The terminal fall speed is also shown (line). 


where Z> mass is the mean mass diameter and p w is the 
density of water. The mean mass diameter is directly 
calculated from observed DSD and is related to the 
slope and shape parameters of the gamma distribution 
as 


D 


mass 


4 + 177 
A 


( 9 ) 


N(D) = N*f l (m ) 


N(D)=NJ 2 (m) 


D 


Z) 


D 


D_ 


exp 


exp 


— (4 + 777 ) 
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— (4 + m) 


Z) 


D 


£> 


and 

( 10 ) 

, ( 11 ) 


where /i (777) and/ 2 (777) are given as 


The normalized intercept parameters can then be cal- 
culated from observed spectra as well. The corresponding 
normalized gamma-fitted distributions are expressed as 




(4 + m) m+1 

r (777 + 1) 


and 


( 12 ) 
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Fig. 9. Mean (squares) and standard deviation (vertical bars) of 
the fall speed measured by PM Tech PD as a function of drop di- 
ameter. The measurements were taken during C3VP and terminal 
fall speed is also given (line). 


/ 2 ( m ) 


6 ( 4 + m ) m+ 4 
256 T(rn + 4) ' 


(13) 


The shape parameter is extracted by minimizing the 
error in rain rates that are calculated from observed 
spectra and from fitted-gamma distribution in Eqs. (10) 
and (11). The shape parameters that are derived from 
Eqs. (10) and (11) differ from each other, and they are 
denoted as nix and m 2 for the rest of the study. It should 
be noted that both Eq. (9) and the rain rate calculated 
from Eqs. (10) and (11) assume the complete gamma 
function where the integral over drop diameter is from 
0 to infinity. In reality, there is a minimum and maximum 


drop size in an observed DSD and the incomplete 
gamma function should be a more realistic fit. The dif- 
ferences in gamma-fitted parameters due to complete 
and incomplete gamma function have been studied in 
the past (Ulbrich and Atlas 1998; Vivekanandan et al. 
2004) and are beyond the scope of this study. 

To determine the measurement accuracy, a subset of 
1-min observations was created from JWD and PD ob- 
servations where at least 100 drops were recorded. At 
least 50 drops were counted for the subset in 2DVD, 
which has twice the sampling cross section. The subset 
had 5349 l-min samples for all DSD and rain parameters 
except for the shape parameter. The range of shape 
parameters may be considered to fall between —2 and 20 
even though higher values are occasionally reported 
(Tokay and Short 1996). Approximately 7% of the data 
was disregarded because it was outside this range. The 
out-of-range values of the shape parameters are partly 
due to the use of complete gamma function and partly 
because the gamma fit is not well suited to the observed 
DSD. 

Four statistics [Eqs. (2)— (5)] are listed in Table 1 for 
three pairs of eight DSD and integral rain parameters. 
The biases show that PD had the highest values and the 
JWD had the lowest values for all parameters. Because 
of the logarithmic values, this may not be clear in N* and 
N w . The percent bias between PD and the other two 
disdrometers exceeded 20%, largely a result of PD’s 
overestimation of drops larger than 2.4 mm in diameter. 
The bias and absolute bias in reflectivity was 0.9 and 
2.0 dBZ between PD and JWD because of the differences 
in midsize and large drops. The differences between the 
disdrometer measurements were also noticeable in the 
mean mass diameter and liquid water content. Conse- 
quently, the intercept and shape parameters of the 
gamma distribution that are based on integral rain and 
DSD parameters showed noticeable differences. While 
the significant differences in small drops were the main 
cause for observed differences in (V*, the differences at 


Table 1. Four different rainfall statistics of the DSD and integral rain parameters between the pairs of disdrometer measurements. 


Bias Absolute bias Bias (%) Absolute bias (%) 



2DVD, 

PD 

2DVD, 

JWD 

PD, 

JWD 

2DVD, 

PD 

2DVD, 

JWD 

PD, 

JWD 

2DVD, 

PD 

2DVD, 

JWD 

PD, 

JWD 

2DVD, 

PD 

2DVD, 

JWD 

PD, 

JWD 

R (mmh 3 ) 

-0.31 

0.17 

0.48 

0.91 

0.73 

1.09 

-7.5 

4.5 

11.9 

21.6 

18.5 

26.6 

Z (dBZ) 

-0.64 

0.25 

0.89 

1.76 

1.38 

2.00 

-2.2 

0.8 

3.0 

6.0 

4.8 

6.9 

^mass (mm) 

-0.09 

0.03 

0.13 

0.13 

0.09 

0.16 

-6.7 

2.7 

9.4 

9.7 

6.6 

11.8 

w (gnT 3 ) 

-0.001 

0.006 

0.008 

0.036 

0.033 

0.044 

-0.7 

3.0 

3.7 

17.5 

16.3 

21.3 

lOlogiV* 

2.39 

-0.90 

-3.29 

2.41 

1.64 

3.34 

10.4 

-3.6 

-14.0 

10.5 

6.7 

14.3 

lOlogW 

1.39 

-0.55 

-1.95 

1.63 

1.13 

2.12 

3.9 

-1.5 

-5.4 

4.6 

3.1 

5.9 

m\ 

-2.51 

0.44 

2.95 

2.81 

1.69 

3.41 

-45.2 

10.7 

55.3 

50.4 

41.4 

63.9 

m 2 

-0.70 

0.85 

1.55 

2.08 

1.77 

2.75 

-13.3 

19.0 

32.2 

39.7 

39.8 

57.2 
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all sizes resulted in significant differences in nil as well. 
The differences in N w and in m 2 were relatively less since 
the differences in the small drop range played a less 
important role for these parameters. It should be noted 
that the differences in drop counts at sizes larger than 
4.5 mm in diameter typically do not contribute to the 
differences in DSD and rain parameters significantly 
because of their infrequent occurrence in disdrometer 
observations. 

The probability and cumulative distributions of R , Z, 
w, and D mass did not show significant differences be- 
tween the disdrometer measurements (Figs. lOa-h), but 
the mean values were noticeably higher in PD and lower 
in the JWD (Table 2). The difference in the mean values 
between the these two disdrometers were 0.5 mm h -1 in 
R , 3 dBZ in Z, 0.008 gm~ 3 in w, and 0.13 mm in D mass . 
The mean values in the 2DVD data were between the 
other two disdrometers, and the median values were 
more similar among the three disdrometers. 

The probability and cumulative distributions of the 
intercept and shape parameters showed significant dif- 
ferences between the disdrometer measurements (Figs. 
lOi-p). The lack of small drops in PD shifted the prob- 
ability and cumulative distributions of Nf and N w to- 
ward smaller values but more so for distributions of Nf. 
The significantly lower mean and median values of Nf 
and N w in PD were reflected in the probability and cu- 
mulative distributions (Figs. 10i,j). The 2DVD also had 
noticeable differences from the JWD in the probability 
and cumulative distributions as a result of the lack of 
small drops in the 2DVD DSDs. The shape parameters 
had noticeably different probability and cumulative 
distributions, but the PD distributions exhibited sub- 
stantially different distributions than the other two dis- 
drometers. The median values of mi and m 2 were 3.2 and 
3.4 in JWD and were very close to the constant value of 
m = 3, which is used in the TRMM precipitation algo- 
rithm (Kozu et al. 2009a). 

To determine the role of the differences in midsize 
and large drops on probability and cumulative distri- 
butions of intercept and shape parameters, these pa- 
rameters were recalculated after eliminating small drops 
from the drop size spectra. The agreement between the 
disdrometers in distributions of Nf and N w was very 
good, but the distributions were quite different than when 
the small drops were included (Figs. lla,b,e,f). The 
agreement in distributions of mi and m 2 was also greatly 
improved when the small drops were eliminated, but 
there were still noticeable differences (Figs. llc,d,g,h). 
The mean and median values of PD were relatively lower 
than the other disdrometers, and all three resulted in 
higher shape parameters and lower intercept parameters 
in the absence of small drops (Table 2). 


9. Conclusions 

This study presented a series of measurement com- 
parisons between JWD, 2DVD, and OTT PARSIVEL. 
The results of the instrument comparisons can be used 
to address questions related to disdrometer application 
and measurement uncertainty. For example, if the re- 
quirement is accurate sampling of the entire drop size 
spectrum, what is the most accurate disdrometer con- 
figuration to employ? The answer obtained from this 
dataset is a combination of the JWD and 2DVD. For the 
study period examined herein, the JWD performed 
relatively well at the small drop end, while the 2DVD 
was able to better distinguish the size of large drops. If 
the requirement is related to the integral parameters of 
w and R, both the JWD and 2DVD can be advocated 
since the small drops do not significantly influence these 
parameters. The contribution of drops larger than 4.5 mm 
in diameter on w and R can also be considered insig- 
nificant since these drops are typically rare in most DSD 
observed with disdrometers. For Z, the very large drops 
can play a significant role especially in the presence of 
low drop counts, which are typically observed at the very 
early stage of convective rainfall in disdrometer mea- 
surement. At the same time, the large drop counts can 
be significant in a weather radar volume of continental 
convection having high values of reflectivity. Thus, the 
significance of large drops on the integral rainfall pa- 
rameters depends upon the measurement scale. 

Concerning fall velocity, the OTT PARSIVEL mea- 
surements followed the parametric representation of the 
GK terminal fall speed of raindrops fairly well, an im- 
provement over previous PM Tech PARSIVEL mea- 
surements that exhibited significantly slower fall velocities 
for drop sizes larger than 1.5 mm in diameter. The im- 
proved fall speed measurements by the OTT PARSIVEL 
seem to be associated with a more sophisticated calibra- 
tion procedure. The 2DVD, on the other hand, recorded 
lower fall velocities than GK-predicted terminal fall 
speeds for drop sizes exceeding diameters of 2.44 mm in 
several events where the mean wind speeds exceeded 
4 ms -1 . The cause of slower fall speeds measured by the 
2DVD is a topic for investigation but beyond the scope of 
this study. 

For parametric forms of the DSD, normalized gamma 
function parameters such as N w can be extracted more 
accurately than those with respect to Nf since the latter 
is more sensitive to the small drop counts. This is im- 
portant because the parametric form of the DSD is of 
great interest for the precipitation retrieval community. 
The integral parameters derived from the parametric 
form of DSD require information on the minimum and 
maximum drop diameters. The minimum drop diameter 
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Fig. 10. Probability and cumulative distributions of (a),(e) R, (b),(f) Z, (c),(g) w, (d),(h) Z) mass , (i),(m) intercept of gamma-fitted 
distribution with respect to N T and D mass , (j),(n) intercept of gamma-fitted distribution with respect to w and D mass (k),(o) shape pa- 
rameter of gamma-fitted distribution with respect to IV*, and (l),(p) shape parameter of gamma-fitted distribution with respect to N„. The 
shape parameter of 3, which is used by the TRMM precipitation radar algorithm, is marked in (k) and (1). 
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Table 2. Mean and median of the distributions of the DSD and integral rain parameters, which are derived from three different 
disdrometer measurements. The values for raindrop distribution parameters that were recalculated after removal of small drops are 
marked in boldface. 




Mean 



Median 


2DVD 

PD 

JWD 

2DVD 

PD 

JWD 

R (mm h~‘) 

4.07 

4.35 

3.86 

2.26 

2.26 

2.23 

Z (dBZ) 

35.3 

37.72 

34.76 

18.76 

18.53 

18.57 

w (gm~ 3 ) 

0.208 

0.210 

0.202 

0.135 

0.128 

0.136 

D mass (mm) 

1.35 

1.44 

1.31 

1.28 

1.36 

1.24 

Nf( m~ 3 mm 1 ) 

390 

181 

419 

239 

145 

305 

N w (m _3 mm _1 ) 

7529 

4355 

8446 

4109 

3115 

4795 

m\ 

4.3 

6.8 

3.9 

3.8 

5.8 

3.2 

m 2 

4.9 

5.6 

4.1 

4.2 

4.6 

3.4 

Nf(m 3 mm 1 ) 

67 

57 

65 

50 

44 

51 

N w (m _3 mm _1 ) 

1770 

1422 

1754 

1307 

1091 

1330 

YYl\ 

7.9 

7.1 

8.6 

7.8 

7.2 

8.7 

m 2 

10.1 

9.2 

10.7 

9.6 

8.7 

10.5 


is mainly linked to the disdrometer’s ability to measure 
the smallest drops. At the current stage, it is considered 
0.2 mm for the 2DVD and 0.3 mm for the JWD. Since 
the 2DVD severely underestimates the number of drops 
in the first size bin, 0.3 mm should be considered the 


minimum drop diameter. On the other end of the size 
spectrum, the maximum drop diameter is related to the 
sampling area of the disdrometer as well as the physics 
of precipitation. Since a radar volume is much greater 
than the measurement area of a single disdrometer, 
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Fig. 11. Probability and cumulative distributions of (a),(e) intercept of gamma-fitted distribution with respect to A r and Z) mass , (b),(f) 
intercept of gamma-fitted distribution with respect to w and Z) mass , (c),(g) shape parameter of gamma-fitted distribution with respect to 
Nj-, and (d),(h) shape parameter of gamma-fitted distribution with respect to N w after eliminating small drops from size spectra. 
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long-term observations with clustered or networked 
2DVDs are needed to determine the maximum drop 
diameter and, subsequently, its impact on the retrieval 
of DSD moments. The GPM ground validation pro- 
gram, which has recently operated five of the third- 
generation 2DVDs in field campaigns and begun an 
extended deployment of them at WFF, will be expected 
to provide information on this matter. 

Underestimation of the number of drops less than 
0.76 mm in diameter and overestimation of the drop 
concentration above 2.4 mm in diameter are the short- 
comings of OTT PARSIVEL. These trends were not 
observed in the PM Tech PARSIVEL. This apparent 
limitation of the OTT PARSIVEL tended to occur 
when mean hourly rain rates exceeded 2.5 mm h _ 1 and 
the total number of drops were greater than 400 in a 
sampled minute but no dependence on environmental 
wind speed was found in this study. Therefore, the ob- 
served measurement bias is attributed to the partially 
inhomogeneous laser beam of the OTT PARSIVEL. 
The manufacturer released an upgraded version of 
the OTT PARSIVEL (PARSIVEL 2 ) where the laser 
beam is expected to be homogeneous. An evaluation of 
PARSIVEL 2 is underway and is expected to perform 
better at all size ranges. 
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